Coexistence of superconductivity and ferromagnetism in two dimensions 
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Ferromagnetism is usually considered to be incompatible with conventional superconductivity, as 
it destroys the singlet correlations responsible for the pairing interaction. Superconductivity and 
ferromagnetism are known to coexist in only a few bulk rare-earth materials. Here we report evi- 
dence for their coexistence in a two-dimensional system: the interface between two bulk insulators, 
LaAlOs (LAO) and SrTiOs (STO), a system that has been studied intensively recently. Magnetore- 
sistance, Hall and electric- field dependence measurements suggest that there are two distinct bands 
of charge carriers that contribute to the interface conductivity. The sensitivity of properties of the 
interface to an electric field make this a fascinating system for the study of the interplay between 
superconductivity and magnetism. 

PACS numbers: 73.20.-r, 73.40.-c, 74.78.Fk, 73.21.-b 



There has been much interest recently in the conduct- 
ing interface that forms between the two band insulators, 
SrTi0 3 (STO) and LaA10 3 (LAO) (Ml- This interface 
shows a rich variety of behavior, including superconduc- 
tivity d, 0, E3, magnetism |jJE EH EH , and electric 
field controlled metal-insulator [2|, [6[ and superconductor- 
insulator transitions 0, [s[ . Two important mechanisms 
have been proposed for the creation of the conducting 
layer at the interface [H, [H, T3-17|: charge transfer from 
the LAO to the Ti +2 ions at the interface (the so-called 
"polar catastrophe" mechanism); and conduction due to 
oxygen vacancies, which can be controlled by growth or 
post-growth annealing conditions. The Ti bands are also 
thought to contribute to the magnetism seen in some 
samples [14j . 

The electronic characteristics are very sensitive to the 
growth conditions: generally, it is found that samples 
grown in an environment with low oxygen partial pres- 
sure Po 2 have more oxygen vacancies and are conse- 
quently more conducting; the conductivity is reduced if 
the samples are grown in high Po 2 , or subjected to a post- 
growth oxygen anneal [si E3 . It is not only the conduc- 
tivity that is sensitive to growth conditions: supercon- 
ductivity is observed in samples grown in intermediate 
Po 2 |S|, and signatures of ferromagnetism are observed 
for samples grown in high Po 2 [HEM, although both phe- 
nomena have not been observed until now in the same 
sample. Here we report measurements on LAO/STO in- 
terface structures where both phases are seen simultane- 
ously at low temperatures. Magnetoresistance and Hall 
measurements indicate that there are at least two bands 
of charge carriers in the system. Earlier theoretical calcu- 
lations [18j have pointed to a ferromagnetic ground state 
of the system due to the multivalent nature of the Ti 
ions at the interface, but the origin of superconductiv- 
ity is still not clear. This system joins only a few other 
bulk materials in which superconductivity and ferromag- 



netism have been observed simultaneously [191-122]. with 



two critical differences: both the superconductivity and 
the magnetism are confined to a two-dimensional inter- 
face, and the electronic properties of this interface can 
be tuned over a wide range by means of an electric field. 
Consequently, it forms a unique system for the investi- 
gation of the interplay between superconductivity and 
magnetism. 

The films in this work had 10 unit cells (uc) of LAO 
grown by pulsed laser deposition at Po 2 = 10 -3 mbar 
on Ti02 terminated (001) STO single crystal substrates 
[23j| . For electrical measurements, a Hall bar geometry 
was patterned by photolithography and etched using ar- 
gon ion milling that removed the LAO layers and a few 
layers of STO. Measurements confirmed that the bare 
etched STO was not conducting from room temperature 
down to millikelvin temperatures. A gate voltage V g ap- 
plied to the back of the substrate was used to modulate 
the conductance of the devices. While the qualitative 
behavior of all the samples was the same, the devices 
showed small sample-to-sample variations that were only 
evident at millikelvin temperatures. The origin of these 
variations is not clear at the moment. We shall concen- 
trate in this paper on data from a single longitudinal 
section and its adjacent Hall configurations that showed 
the sharpest superconducting transitions. 

The normal state and superconducting characteristics 
and their dependence on V g and temperature T are sim- 
ilar to that seen by other groups [jj 0, [s[ (Fig. 1(a)). 
The current- voltage curves also show a characteristic su- 
perconducting signature with a critical current I c that 
vanishes when V g < —20V (Fig. 1(b)). However, I c and 
the transition temperature T c do not increase monotoni- 
cally with increasing V g , but show maxima at V g ~ 80 V 
(Fig. 1(c)). 

Figure 2(a) shows the superconducting transition at 
V g = 80 V at a few different values of magnetic field H 
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FIG. 1: (a) Superconducting transition at different gate volt- 
ages V g . (b) Current voltage characteristics of the sample at 
15 mK. (c) Transition temperature T c (defined as the mid- 
point of the resistive transition) and critical current (defined 
as the current at which the sample switches to the resistive 
state on ramping the current up from zero) as a function of 
V g . Both measures indicate that the maximal superconduct- 
ing properties are obtained for V g ~ 80 V. Error bars indicate 
the difference in T c measured between cooling and warming 
traces. 

applied perpendicular to the film plane. Defining T c as 
the temperature corresponding to half the normal state 
resistance, a plot of T C {H) is shown as an inset to the 
figure. For two dimensional superconductors in this field 
orientation, T C (H) should be linear, its slope proportional 
to the Ginzburg-Landau coherence length £o [24] . A lin- 
ear fit gives £o = 64 nm, close to that reported previously 

SEo]. 

However, this method of determining T C (H) misses 
some of the more novel and interesting behavior of this 
system. To show this behavior, we map T C (H) continu- 
ously by controlling the sample temperature so that the 
sample resistance Rs remains fixed as we ramp the mag- 
netic field. Figure 2(b) shows T C (H) for three different 
resistance bias points at V g = 80 V. The most striking 
aspect of the data is that the behavior is hysteretic as 
a function of H. Consider the curve corresponding to 
the bias point Rs = 208 fi, slightly lower than the mid- 
point of the transition. Increasing the magnetic field from 
negative values towards H = 0, T c shows a smooth in- 
crease that is almost linear. The slope of this linear curve 
gives £o = 71 nm, close to the value obtained from Fig. 
2(a). On ramping H beyond zero, however, the behavior 
becomes non-monotonic: in particular, there is a local 
minimum at H ~ 7 mT. Increasing H further, T C {H) 
becomes monotonic again. Reversing the ramp direction 
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FIG. 2: (a) Superconducting transition for V g = 80 V, at 
a few different magnetic fields. Inset: T c vs. H obtained 
from these data. T c is defined by the midpoint of the resistive 
transition. The line represents a linear fit to the data points, 
with a slope of 1.6 mK/mT, giving a zero temperature coher- 
ence length of = 64 nm. (b) Phase diagram, T c vs. if, 
for a gate voltage of V g = 80 V, where the superconducting 
properties are maximal. The three curves represent different 
resistance bias points along the superconducting transition, 
with the normal state resistance per square being 704 Q. The 
arrows mark the direction of the magnetic field ramp, (c) 
Phase diagrams at a bias point of Rs = 83 Q, at the foot of 
the resistive transition, at three different gate voltages. Data 
for V g = 40 V are shifted by 10 mK and data for V g = 80 V 
by 20 mK for clarity. 

results in a mirror image of the first trace, giving rise 
to a characteristic "butterfly" curve. Similar behavior 
is seen at other gate voltages for which the sample goes 
superconducting (Fig. 2(c)). 

Hysteretic behavior and butterfly curves in the mag- 
netic field dependence of electrical characteristics are sig- 
natures of underlying ferromagnetic order in a sample. 
Such behavior has already been reported in LAO/STO 
interface samples: Brinkman et al. [4] demonstrated that 
hysteretic behavior is observed at low temperatures in 
26 uc thick LAO films grown under high Po 2 \ more re- 
cently, Ariando et al. [13] were able to observe hysteretic 
magnetization curves coexisting with paramagnetic and 
diamagnetic behavior in samples also grown in high Po 2 
that persisted to room temperature. For low Po 2 during 
growth, Huijben et al. [8[ note that a low sheet resis- 
tance and metallic behavior is observed; for intermediate 
Po 2 , the samples go superconducting; while for high Po 2 
magnetic behavior is seen. These experiments suggest 
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to us that the magnetism is associated with the polar 
catastrophe mechanism, while the superconductivity is 
associated with the presence of oxygen vacancies. How- 
ever, in all previous experiments, the magnetic and su- 
perconducting regimes were quite distinct. In contrast, 
our samples, which are grown under high Po 2 and have 
Rs consistent with the superconducting samples mea- 
sured by others [8[, show a coexistence of superconduc- 
tivity and ferromagnetism. The electronic properties of 
LAO/STO interfaces are extremely sensitive to growth 
conditions: samples grown by different groups with nom- 
inally identical growth conditions show some variations 
in electronic properties. Our samples appear to be in a 
growth regime where both phenomena coexist. 

Caviglia et al. [l2[ have reported weak (anti- 
)localization magnetoresistance (MR) in LAO/STO sam- 
ples. At first sight, the data from our samples appear 
to be very similar to their data. Figure 3(a) shows the 
MR of the sample at a few different temperatures at 
V g = —100 V, not in the superconducting regime. The 
magnitude of the resistance change and the sharpness of 
the resistance dip near zero field increase with decreasing 
temperature, consistent with a phase coherence length 
that increases with decreasing temperature. However, a 
closer look at the low field MR, shown in Fig. 3(b), re- 
veals some significant differences from [12| . In addition to 
the nonmonotonic MR over a large field scale (Fig. 3(a)), 
the samples show an additional resistance dip near zero 
field whose magnitude increases and width decreases with 
decreasing temperature. This behavior, which is char- 
acteristic of weak localization (WL) [25j, suggests that 
there are two independent carrier gases that contribute 
to the conductance of the device in parallel, each with its 
own WL contribution. The MR also shows a hysteretic 
"butterfly" pattern similar to that seen in T C (H) (Fig. 
2(b,c)), indicating that local magnetic fields arising from 
magnetic order also modulate quantum interference in 
the carrier gases. 

To demonstrate the feasibility of this picture, we 
have simulated the WL contribution of two parallel two- 
dimensional gases in the presence of an external magnetic 
field and a magnetic field arising from a hysteretic intrin- 
sic magnetization. To calculate the WL contribution in 
the presence of spin-orbit scattering, we use the formal- 
ism of Santhanam et al (26|. The resulting curve is 
shown in Fig. 3(c). The simulation qualitatively repro- 
duces the experimental features: a non-monotonic MR 
over a large field scale, and a smaller non-monotonic MR 
over a smaller field scale, with hysteresis due to intrinsic 
magnetic order. In the simulation, the larger contribu- 
tion (88%) is due to a carrier gas that has short phase 
coherence (L^) and spin-orbit scattering (L so ) lengths 
(L4 = 0.17 /im, L S o — 0.03 /im), while the remaining 
contribution is due to a carrier gas with longer and 
L so = 4.5 //m, L so = 0.45 /mi). Experimentally, we 
observe an increase in the amplitude of the low-field MR 
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FIG. 3: (a) Magnetoresistance (MR) at V g = -100 V at four 
different temperatures, (b) Expansion of the low field data 
of (a) for the three highest temperatures, (c) Simulation of 
the longitudinal MR assuming weak localization (WL) con- 
tributions from two parallel conducting channels. The inset 
shows an expanded view of the zero field region, showing the 
hysteresis and the non-monotonic dependence seen in the ex- 
perimental data, (d) Hall resistance at V g = — 100V at three 
different temperatures. Data for 50 mK and 375 mK have 
been offset by + 5 and -5 O respectively, (e) Hall coefficient 
and zero field longitudinal resistance at different gate volt- 
ages at 270 mK. The longitudinal resistance decreases with 
increasing gate voltage, but the Hall coefficient, nominally 

inversely proportional to the electron density, increases. 

as the sample approaches the superconducting transition 

(by changing V g , for example), which is probably due 

to the contribution of Maki-Thompson superconducting 

fluctuations |27j whose MR field scale is determined by 

[26(| , and which is expected to increase exponentially 

as T T c . 

Further evidence for two parallel conduction channels 
can be seen in the Hall resistance data. Figure 3(d) shows 
the Hall resistance RhclII of the sample at V g = —100 V 
at three different temperatures. The sign of the slope of 
RHaii(H), Rh, corresponds to negatively charged carri- 
ers with a density n ~ 4.4 x 10 13 /cm 2 , assuming a single 
carrier band. If there were only one electron gas at the in- 
terface, increasing V g should increase n, decreasing both 
Rs and Rh- As shown in Fig. 3(e), increasing V g does 
indeed decrease Rs, but increases Rh- This behavior is 
only possible if there are at least two types of carriers in- 
volved in electrical transport, with different dependences 
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FIG. 4: (a) Magnetoresistance (top curve) and Hall resistance 
(bottom curve) at V g = 20 V at 50 mK. (b) Same data as 
in (a), for one direction of the magnetic field sweep, except 
the background linear component is subtracted from the Hall 
resistance. 

of densities and mobilities on gate voltage. Evidence for 
multiple charge carriers has also been observed by other 
groups [II [l7j]. 

For conventional itinerant ferromagnets, it is well 
known that a finite magnetization should manifest it- 
self as a contribution to RHaii through the anomalous 
Hall effect [28| . Figure 4(a) shows the longitudinal and 
Hall MR at V g = 20 V and 50 mK, where R s =0 when 
H=0 (Fig. 1(a)). Figure 4(b) shows the same data with 
the field sweep in one direction. Here a linear contri- 
bution determined by fitting RHaii(H) at high magnetic 
fields has been subtracted. It can be seen that some of 
the structure in Rs(H) also has corresponding signatures 
in RHaii(H) (this structure is not due to misalignment 
of the Hall contacts, which is very small). The struc- 
ture in Rs(H) and R Ha ii(H) for H > (for this field 
sweep direction) is also seen at other gate voltages, fur- 
ther from the superconducting transition. We believe 
that this structure in the Hall resistance is due to an 
anomalous Hall effect arising from the interaction of the 
charge carriers with the magnetic moments at the inter- 
face, but the structure seen in Fig. 4(b) indicates that 
this interaction is more complicated than that in simple 
itinerant ferromagnets, and warrants further investiga- 
tion. The features in Rho,ii(H) change with V g . The 
sharp dip in both the longitudinal and the Hall resis- 
tance for H < appears only close to the superconduct- 
ing transition, and may be associated with vortex flow in 
the superconductor [29j. 

In summary, the interface between LAO and STO 
shows a rich variety of behavior, including interacting fer- 
romagnetism and superconductivity. Our measurements 
indicate that there are two different types of charge car- 



riers in the system. The ability to tune the properties 
of the system by means of a gate voltage makes this a 
fascinating system for studying competing cooperative 
phenomena in two dimensions. 
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